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1 INTRODUCTION 



In a recent paper llTownsendll2005l. hereafter T05), I inves- 
tigated the influence that the rotation-originated Coriolis 
force exerts over the self-excited oscillations of slowly pul- 
sating B-type (SPB) stars. To maintain a link with his- 
torical theoretical studies of SPB stars in which the ef- 
fects of rotation were overloo ked (e.g., iGautschv fc Saiol 
119931: iDziembowski et ai]|l993ft . the analysis was restricted 
to those modes that transform continuously into ordinary g 
modes as the rotation rate is reduced toward zero. Following 
the classificati on scheme for low-frequency waves in rotat- 
ing fluids fe.g.. lGilll ll982-"l. these comprise the prograde and 
retrograde Poincare modes, the prograde mixed modes, and 
the Kelvin modes. 

The present paper extends the T05 study to one of the 
types of pulsation mode found only in rotating stars: the 
retrograde mixed modes, also referred to in the l iterature 
as mixed Rossby-Poincare modes (e.g , ILouI |200cT) . mixed 
Ross by-gravity modes (e.g., lGilllll982i) . and Yanai modes 
fe.g.. lTownsendll2003aft . After providing a brief background 
to these modes in Section [5] I use an approximation-based 
method (Sec. [3J, very similar to that developed in T05, to 
investigate their stability in a range of B-type stellar mod- 
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ABSTRACT 

I examine the stability of retrograde mixed modes in rotating B-type stars. These 
modes can be regarded as a hybrid between the Rossby modes that arise from con- 
servation of vorticity, and the Poincare modes that are gravity waves modified by the 
Coriolis force. Using a non-adiabatic pulsation code based around the traditional ap- 
proximation, I find that the modes are unstable in mid- to late-B type stars, due to 
the same iron-bump opacity mechanism usually associated with SPB and (5 Cep stars. 
At one half of the critical rotation rate, the instability for m = 1 . . . 4 modes spans the 
spectral types B4 to AO. Inertial-frame periods of the unstable modes range from 100 
days down to a fraction of a day, while normalized growth rates can reach in excess of 
10~ 5 . 

I discuss the relevance of these findings to SPB and pulsating Be stars, and to the 
putative Maia class of variable star. I also outline some of the questions raised by this 
discovery of a wholly-new class of pulsational instability in early-type stars. 

Key words: instabilities - stars: oscillation - stars: rotation - stars: variables: other 
- stars: early-type - stars: emission-line, Be 



els (Sec. . The principal finding (Sec. |SJ is that retrograde 
mixed modes are unstable in mid to late B-type stars, due 
to the same iron-bump « mechanism responsible for the pul- 
sation of SPB and /3 Cep stars. I discuss and summarize the 
significance of this result in Section [5] 



2 BACKGROUND 

In order to understand the nature of retrograde mixed 
modes, and why they might be susceptible to K-mechanism 
excitation, it is helpful to review their relationship to the 
differing classes of low-frequency mode present in rotat- 
ing stars. These classes comprise the Poincare, Kelvin, and 
Rossby modes on the one hand, and the mixed modes (both 
prograde and retrograde) on the other. 

The Poincare modes rely on a combination of buoy- 
ancy and the Coriolis force to restore displaced fluid ele- 
ments to their equilibrium position. They reduce to ordi- 
nary gravity modes in the non-rotating limit, having an an- 
gular dependence described by the spherical harmonics Yf m 
with harmonic degree £ and azimuthal order m satisfying 
— t + 2 ^ m ^ I . In light of this correspondence, it is com- 
mon to refer to the Poincare modes sim ply as g modes (e.g., 
iLee fc Saidll987l: iBildsten et alll 1996ft . 

However, such a designation groups the Poincare modes 
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with the Kelvin modes, since the latter also reduce to g 
modes, having m = —£, in the non-rotating limit. These 
Kelvin modes, which always propagate in the prograde az- 
imuthal direction 1 , have rather different properties than the 
Poincare modes, making it desirable to maintain the distinc- 
tion between the two. Most significantly, Kelvin modes ex- 
hibit gcostrophic balance, where the polar component of the 
Coriolis force remains in balance wi th the forc e arising from 
polar pressure gradients (see, e.g.. lGililll982t) . As a result, 
the horizontal propagation of these modes is dispersion-free, 
distinguishing them from the dispersive Poincare modes. 

The Rossby mode s (or r modes; see 
iPapaloizou fc Pringld Il97ct) are unlike the Poincare 
and Kelvin modes, inasmuch as they do not require 
buoyancy to restore displaced fluid elements. Instead, they 
rely on the conservation of the radial component of total 
vorticity, operating on fluid displacements acros s the curved 
level surfaces of the star's stratification. [See ISaid (Il982f) 
for a particularly lucid description of this process]. Rossby 
modes always propagate in the retrograde direction, and in 
the limit of slow rotation follow the first-order dispersion 
relation 



2mQ. 
£(£ + 1) 



(1) 



independent of the internal stellar structure. Here, and 
throughout the present work, the symbols have the same 
meaning as in T05; in particular, a is the pulsation angular 
frequency in the co-rotating frame of reference, and Q is the 
rotation angular frequency. In the same slow-rotation limit, 
the horizontal velocity fields generated by r modes approx- 
imate those produced by the trivial toroid al modes found 
in non-rotating stars (e.g.. lUnno et al1 [l989) . with an angu- 
lar dependence described by the curl of a radial vector with 
length Y e m . 

Mixed modes can be regarded as a hybrid between 
Rossby and Poincare modes. In the dispersio n diagram 
for low- frequency modes in a rotating star (see lLovJl200Ct 
Fig. 1), they mark the dividing line between the Poincare 
modes at higher a, and the Kelvin and Rossby modes at 
lower a. At slow rotation rates, the prograde mixed modes 
behave like Poincare modes, and are usually classified as 
the £ — —rn + 1 g modes. In contrast, the retrograde mixed 
modes have a Rossby-like character, and are classified as r 
modes described by the dispersion relation Q with I — m. 
As the rotation rate increases, however, the prograde and 
retrograde mixed modes for each value of m approach one 
another in frequency, and exhibit properties that lie some- 
where between pure Rossby and Poincare modes. 

There is a well-established body of literature on mixed 
modes in the field of terrestrial atmospheric and oceanic 
geophysics, exte nding back to their d e tectio n in the Earth's 
atmosphere by lYanai fc Maruv amal (Il966l) . coupled with 
the development~bv HM!feundll^)6ai of the 'equatorial 
/3-plane' approximation for analytical treatment of low- 
frequency waves in rotating fluids. Literature discussing 
mixed modes in a stellar astrophysica l context is, how- 
ever, comparatively scarce. ILee fc Said <1987t Il997f) and 



iDziembowski fc Kosovichevl 1119871) have both noted the 
presence of r-mode solutions to the governing rotating-star 
pulsation equations that acquire the character of g modes in 
the limit of large spin parameter v = 2fi/<r. Yet, a full ap- 
preciation of the nature of these mixed-mode solutions did 
not emerge until the /3-plane approximation was appli ed to 
the st ellar pulsation problem , in the recent studies by 
l|2000|) and lTownsendl i2003aTl . 

From a standpoint of stellar stability, the retrograde 
mixed modes are particularly interesting on account of their 
relati onship to r modes. Being almost solenoidal (e.g.. ISaicI 
1982), r modes produce little compression or rarefaction of 
the stellar plasma, and are therefore difficult to excite via 
the classical 'Gre ek letter' mechanisms (7, 8, e and k; see 
lUnno et alJll989h which involve a thermodynamic Carnot 
cycle. This explains why incidences of r-mode instability due 
to these mechanisms are rare 2 . Indeed, even in those situa- 
tions where classical instability has been found, the strength 
of the excitation - as measured by its linear growth rate - 
has been much smaller than is typical for or dinary p or g 
modes (see, e.g.. iBerthomieu fc ProvosHll983l) . 

However, in the case of retrograde mixed modes the 
non-solenoidal Poincare character acquired with more-rapid 
rotation causes them to generate appreciable compressions 
and rarefactions. If rotation is sufficient, therefore, these 
modes stand a good chance of becoming (significantly) un- 
stable to one of the classical excitation mechanisms. This 
recognition is the foundation for the present paper. 



3 METHOD 

The method used to investigate the stability of retrograde 
mixed modes is based on the so-called 'traditional approx- 
imation', and is very similar to the approach described in 
detail in T05. In particular, I solve an identical set of cou- 
pled differential equations (7-10, ibid) and boundary con- 
ditions (15-16, ibid) as an eigenproblem, using the non- 
adiabatic BOOJUM code. The resulting dimensionless eigen- 
frequencies 



GM, 



(2) 



1 I adopt the usual convention that negative m corresponds to 
prograde propagation, and positive m to retrograde propagation, 
as viewed from the co-rotating frame. 



define, via their imaginary parts aJi, whether a mixed mode 
is linearly stable (ui > 0) or unstable (u>i < 0). The most im- 
portant departure from T05 relates to the calculation of the 
effective harmonic degree I appearing in the pulsation equa- 
tions. In T05, I is obtained (eqn. 6, ibid) from the eigenvalue 
A™ of Laplace's tidal equations that transforms continuously 
into £(£ + 1) in the limit of no rotation. 

In the present analysis, I instead evaluate I from the 
family of eigenvalues A that are appropriate to retrograde 
mixed mod es. Within the alte rnat ive indexing schem es in- 
troduced bv lLee fc Said" Jl997h and lTownsendlll2003al) . they 
are labeled as the k — — 1 or s = solution branches of the 



2 It should be noted that there exist other, non-classical mecha- 
nisms for the excitation of r modes (e.g., driving by gravitational 
waves; see Andcrsson 1998), but these are generally unimportant 
in the non-relativistic objects considered here. 
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tidal equations. Only minor revisions to BOOJUM were re- 
quired to permit consideration of these branches, consisting 
primarily of alterations to the matrix routines employed to 
calculate A. However, initial usage of the updated code re- 
vealed very poor convergence toward eigenfrequencies, with 
thousands of relaxation iterations required to reach the de- 
sired solution tolerance. After some investigation, the dif- 
ficulty was traced to the adoption in T05 of the steady- 
wave approximation. This approximation, which entails us- 
ing only the real part of the frequency when evaluating u, 
permits easy solution of the tidal equations since all vari- 
ables remain real. However, it has the undesirable side effect 
that the implicit dependence of A on ui becomes non-analytic 
(i.e., A does not satisfy the Cauchy-Riemann equations). In- 
variably, this tends to hamper the convergence of BOOJUM, 
especially when considering mixed modes. 

This difficulty is addressed within BOOJUM via a new 
two-stage procedure for evaluation of A. First, Laplace's tidal 
equations are solved as before using only the real part ui r of 
the frequency. The solutions are then taken a s the start- 
ing p oint in an inverse iteration algorithm (e.g.. |Press et all 
1992) based on the fully-complex frequency lj. The algo- 
rithm is allowed to run until the relative change in A reaches 
machine precision, which typically takes five or six itera- 
tions. With this modification, A behaves in an analytic fash- 
ion, and the convergence issues are wholly resolved. 

Further improvements to BOOJUM concern the discrimi- 
nant root finding process. Although the dimensionless eigen- 
frequencies of the retrograde mixed modes remain defined by 
the roots 



£>(«) = 



(3) 



of the discriminant function (see T05, Sec. 2.3), I found that 
initial isolation of roots was greatly facilitated by adopting 
the quantity 



Slj = 



U)Q — UI 

n 3 



(4) 



as the independent va riable , inste ad of u> itself. This quan- 
tity, also favoured bv ISaiol I0~982), furnishes a normalized 
measure of departures from the first-order dispersion rela- 
tion Q. The dimensionless rotation frequency appearing in 
the denominator is defined as 



Q = Q 



m 

GM, 



while the r-mode limit for mixed modes, 



Wo 



20. 
m + 1' 



(5) 



(6) 



comes from setting £ — m in the dispersion relation Q. 
Since deviations from sphericity arising due to the centrifu- 
gal force are neglected in the presen t analysis, u> is always 
smaller than ujo (see, e.g.. ISaiall982h . and thus the 'depar- 
ture measure' b~u is positive for all modes considered. 

The algorithm used by BOOJUM to solve the character- 
istic equation ^ has also been alte red. I found t hat the 
complex secant approac h favou r ed by lCastorl jl97ll) consis- 
tently outperforms the iMullerl l)l956h algorithm originally 
implemented, even though the latter method has a higher 
(quadratic) order of convergence toward roots. 
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Figure 1. The theoretical HR diagram for the stellar models 
introduced in Section Q] these are plotted as points in the effec- 
tive temperature (T c g) versus stellar luminosity (L*) plane. The 
dashed lines running diagonally from top left to bottom right 
indicate the ZAMS and TAMS main sequence limits, while the 
circle is centred on the position of the 53 Per model analyzed 
in Section 15. II Full evolutionary tracks for three selected initial 
masses (labeled at the ZAMS in solar units) are also shown in the 
diagram, drawn as solid lines. 



4 STELLAR MODELS 

I use the Warsaw-New Jersey evolution code to calculate 
tracks of stellar models for a metalicity Z = 0.02. These 
tracks, of which there are 91 together comprising 4268 in- 
dividual models, sample the initial mass range M* = 2.3- 
5.3 Mq at a resolution 0.05 Mq and the range M» = 5.3- 
8.3 Mq at a resolution 0.1 Mq; each extends from zero- 
age main sequence (ZAMS) to terminal-age main sequence 
(TAMS). This coverage of the theoretical Hertzsprung- 
Russel (HR) diagram, plotted in Fig. is shifted toward 
rather lower masses than considered in T05, primarily to 
ensure - based on initial exploratory calculations - that 
complete coverage of the mixed-mode instability strips is 
secured (see Sec. 15.21 . However, apart from the differences 
in initial masses, the parameters and input data I adopt for 
the stellar model calculations are identical to those discussed 
in Section 3 of T05. 



5 STABILITY CALCULATIONS 
5.1 53 Per Model 

I first apply the updated BOOJUM code to investigate the sta- 
bility of to = 1 mixed modes in the same 53 Per-like model 
introduced in T05. A range of rotation angular frequencies, 
from the negligibly rotating case up to the half-critical rate 
Q/f2 c = 0.5, is considered. This upper bound is loosely dic- 
tated by the limitations of the approximations described in 
Sectional (and see also T05), but also happens to be close to 
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Figure 2. Departure measures for m = 1 mixed modes of the 
53 Per model, plotted as a function of rotation frequency. Each 
curve corresponds to a particular radial order n; the line weight 
is used to indicate whether modes are stable (thin) or unstable 
(thick). 



the observed ceiling on SPB-star rotati on, Q/fl c ~ 0.46 3 , es - 
tablished by the B5 pulsator HD 1976 jMathias et al.l200ll) . 

The results from these stability calculations are plotted 
in Fig. [5] which traces the Q dependence of the eigenfre- 
quency departure measure 5ui introduced in eqn. 0, for 
modes of differing radial orders 4 n. At slowest rotation, 
the flatness of the Slu curves indicates that the departure 
measure f or each mode is independent of Q, confirming the 
finding by P apaloizou fc Pringld 1(1978) that deviations from 
the first-order dispersion relation Q vary as 0(fi 3 ). Toward 
more-rapid rotation, higher-order terms become significant, 
with the result that the curves shown in Fig. |5] dip down- 
ward. 

Of particular interest in the figure is the fact that a 
group of modes, having radial orders n = 13-25, become 
unstable above a rotation rate f2/f2 c « 0.3. The nature of 
the instability is revealed in Fig. [3] which plots both the 
cumulative work W and differential work dW/dx for the 
n — 20 mixed mode, at a rotation rate Q/Q c — 0.25 when 
the mode is stable, and at the more-rapid rate f2/f2 c = 0.5 
when the mode has subsequently been destabilized. The 
sharp peak in dW/dx at log T « 5.3 reveals that the mode is 
strongly driven by the same iron-bump opacity mechanism 
consid ered responsible for the p-mode instability of /3 Cep 
stars dCox et all Il992t IDziembowski fc Pamiatnvkbl 119931) 



3 In fact, this quoted value should be taken as a lower estimate, 
since the inclination of the star remains uncertain. 

4 As enumerated by the number of nodes in the radial displace- 
ment eigenfunction 



and the g-mode instability of SPB stars iGautschv fc Saiol 
Il993t IDziembowski et aTlll993l) . 

As elucidated by the latter-most authors (and see also 
iPamvatnvk h 1999), two key requirements must be satisfied 
for the k mechanism to be effective. First, the relative La- 
grangian pressure perturbation Sp/p should be large and 
vary slowly with radius within the driving zone. This con- 
straint on the shape of the eigenfunction establishes a lower 
limit n ^ 13 on the radial order of unstable mixed modes 
in the 53 Per model, surprising similar to the corresponding 
limit found in T05 for unstable g modes. In fact, the coin- 
cidence between the two stems from the fact that the same 
set of pulsation equations governs both classes of mode. 

The second requirement is that the pulsation period 
should not be significantly greater than the local thermal 
timescale within the excitation zone. If instead the period 
is too long, radiative diffusion between neighbouring fluid 
elements can lead to strong damping in the inner envelope, 
that cancels out any driving due to the k, mechanism. This 
inhibition effect is illustrated in Fig.|3] The strength of the 
K-mechanism driving is the same in both panels, as can be 
seen from the similar peaks dW/dx « 0.04 Ew in the differ- 
ential work at logT ~ 5.3. However, the radiative damping 
at logT ~ 5.7 is more pronounced in the left-hand panel, 
owing to the longer period of the mode - P = 2.76 d, as 
compared to P — 1.63 d for the right-hand panel. This en- 
hanced damping is sufficient to stabilize the mode. 

All mixed modes in Fig.|5|having f2/S7 c < 0.3 or Su > 4 
posses periods longer than the thermal timescale in the driv- 
ing zone. Hence, these modes are stabilized by the same ra- 
diative damping described above. In the following section, I 
explore how radiative damping, acting in tandem with the 
restriction on the 5p/p eigenfunction morphology, shapes 
which regions of the HR diagram exhibit mixed-mode in- 
stability via the K-mechanism. The remainder of the present 
section is focused on a discussion of the unstable n = 1 
mixed mode seen in Fig. [5] as the thick line having 8u ~ 0. 

Inspection of the work function for this mode indicates 
that weak driving in the core, at log T « 7.4, is responsible 
for the instability. Although this may be a result of genuine 
e-mechanism excitation, the normalized growth rate of the 
mode, 

77 = -Ui/u) T , (7) 

never exceeds 3.3 x 10~ 13 at any rotation rate considered 
(compare with the n = 20 mode, which reaches rj = 7.3 x 
10 -6 at Q/Q, c = 0.5). This value differs so little from zero 
that it would be more accurate to describe the n = 1 mixed 
mode as neutrally stable. Such a stance is lent support by the 
fact that the mode sits very close to the pure r-mode limit 
uj = ujo, and - being almost solenoidal - can be expected to 
be neither excited nor damped by classical mechanisms. 



5.2 Instability Strips 

The scope of the investigation is now broadened, to encom- 
pass the entire set of stellar models introduced in Section [1] 
I use BOOJUM to investigate the stability of these models 
toward m = 1 ... 4 retrograde mixed modes. Fig. 0] plots 
the instability strips in the HR diagram associated with 
these modes, at two distinct rotation rates: Q/Q c — 0.5 
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Figure 3. The cumulative work (thin) and differential work (thick) for the m = 1, n = 20 mixed mode of the 53 Per model, plotted as 
a function of temperature at rotation rates f2/Q c = 0.25 (left panel) and 0.5 (right panel). The work functions are expressed in scaled 
units of the total energy of pulsation E\y ■ 



approximates the observed upper limit on SPB-star rota- 
tion (see Sec. 15.2^ . while fi/fi c = 0.25 is representative of 
more-moderately rotating stars. Formally, the condition of 
instability is defined by the requirement that a mode has a 
positive growth rate. However, it is appropriate to exclude 
those modes that, similar to the n = 1 mode considered in 
the preceding section, have such small values of 77 that they 
are more properly regarded as neutrally stable. Therefore, 
the instability strips shown in Fig. [I] are constructed to en- 
close all stellar models that possess one or more modes with 
a growth rate in excess of the threshold value rj = 10 -9 . 

Also shown in Fig. 2] are the inferred positions of se- 
lected SPB and pulsating Be stars, and the position of the 
Maia candidate star HD 208727, which is discussed further 
in Section HJ T he SPB data are comp il ed directly f r om th e 
tabulations by I Waelkens et alJ (I1998T). lAerts et alJ ^199Sf> . 
iMathias et"atl (boOll) and iDe Cat fc Aertsl j20o3T : in the 
many cases where multiple (Terr , L*) values exist for a given 
star, the most recently published ones are adopted. The pul- 
sating Be-star data are obtained via a three-stage process. 
First, the list of line -profile variable Be stars observed by 
iRivinius et"al"l (120031) is merged with that of the short-period 
photometric variables st udied bvlPercv et al.1 l2002tl2004> . 
Then, the calibrations bv lChauville et al.1 J2001T ) are used to 
determine the effective temperature and surface gravity of 
each star in the merged list; if values are not available for 
a given star, it is dropped from the list. Finally, I define an 
empirical relation L* = L*(T c h, log g) via a quintic surface 
fit to the luminosity data of the stellar models introduced in 
Section^] This relation is used to map the (T c g, log g) values 
in the list into corresponding (Tea, L*) values, suitable for 
plotting in the HR diagram. 

Figure 2] demonstrates that the mixed-mode instabil- 



ity found in the 53 Per model (Sec. 15.11 extends to many 
of the other stellar models considered. At a rotation rate 
n/fic = 0.25, the instability for the m — 1 ... 4 modes spans 
the effective temperature range logT c ff = 3.96-4.17. At the 
more-rapid rate fi/n c = 0.5, the instability broadens and 
shifts to toward higher effective temperatures and luminosi- 
ties, with a range log Teff = 3.99-4.23 corres ponding loosely 
to spectral types B4 to AO feohm-Vitenselfr98ll) . At fixed 
rotation rate, incrementing the azimuthal order displaces the 
instability in the opposite direction, toward lower effective 
temperatures and luminosities. 

The sensitivity of the instability toward both rotation 
rate and azimuthal order comes from the period matching 
requirement discussed in the preceding section. Returning 
once again to the first-order dispersion relation Q, it is clear 
that any increase in fl (m) will act to increase (decrease) the 
frequency of a mixed mode, and hence to shorten (lengthen) 
its period. To ensure that this period remains commensurate 
with the local thermal timescale in the driving zone, the in- 
stability therefore shifts to hotter (cooler) stars, in which the 
opacity bump is situated shallower (deeper) in the envelope. 
A very similar line reasoning was employed in T05 to un- 
derstand the behaviour of the instability strips for Poincare, 
Kelvin and prograde mixed modes. 

Examining the correspondence between the mixed- 
mode instability strips and the positions of the selected pul- 
sating stars plotted in Fig. 21 it is clear that there is not a 
great deal of agreement between the two. This discrepancy 
is underlined by Table^ which summarizes the proportion 
of stars enclosed by the instability strips at each azimuthal 
order and rotation rate considered. At best, the strips en- 
compass only two thirds of the SPB stars - a reflection of 
the stars' bias toward the late end of the B spectral type, 
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Figure 4. Instability strips in the HR diagram, for m = 1...4 mixed modes of the stellar models introduced in Section |I] For each 
value of the azimuthal order, the extent of the re-mechanism instability at the rotation rates Q/Q c = 0.25 (dotted) and 0.5 (dashed) is 
indicated using thick lines. The small open triangles, filled squares, and open diamond show the inferred positions of selected SPB stars, 
pulsating Be stars, and HD 208727, respectively (see text). 



whereas SPB stars tend to be concentrated more around 
the mid-B types. In fact, the SPB stars coincide far more 
closely with the theoretical instability strips pr edicted for k- 
mech anism excitation of ordinary g modes (see IPamvatnvkbl 
Il999l. his figs. 3 & 4). Therefore, it seems unlikely that the 
SPB phenomenon as a whole can be explained better by 
mixed modes than by the widely-accepted g-mode interpre- 
tation. 

The situation is superficially similar for the pulsating 



Be stars shown in Fig. [1] in no case do more than 7 percent 
of these stars fall within the mixed-mode instability strips. 
In light of the concentratio n of Be stars toward the early-B 
spectral types JPorterll996l) . and the conspicuous abs ence of 
any pu ls ating Be stars later than sp ectral type B6 (seei Baadc 
1989a|3 IPorter fc Riviniusl I2003L and references therein), 
this result is hardly surprising. However, in contrast to SPB 
stars, Be stars are known to rotate at a significant frac- 
tion of their critical rate ((14/14) « 0.83, corresponding to 
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Figure 5. Maximal growth rates in the HR diagram, for m = 1 . ..4 mixed modes of the stellar models introduced in Section |I] For 
each value of the azimuthal order, the largest rj at the rotation rate Q/Q c = 0.5 is indicated via the gray level. The distribution of these 
levels is logarithmic, and extends from 10 — 9 (white) to 10 — 4,5 (black). 



Table 1. Proportions, expressed as percentages, of the SPB (Be) 
stars that fall within the mixed-mode instability strips plotted in 
Fig-E] f° r each azimuthal order and rotation rate considered. 

m 

n/n c 1234 

0.25 41 (3) 13 (0) 3 (0) 1 (0) 
0.50 67 (7) 48 (3) 41 (3) 29 (3) 



(f2/fi c ) « 0.96; see IChauville et alJ l200lT) . Although near- 
critical rotation lies beyond the reach of the BOOJUM code, 
the general trends shown in Fig. [I] suggest that as fi/fl — > 1 
the mixed-mode instability strips should shift toward the 
early-B spectral types, encompassing a much greater pro- 
portion of pulsating Be stars. Therefore, it remains plau- 
sible that the variability of these stars can be attributed 
to mixed modes, rather than the Poincare modes that are 
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Figure 6. Instability regions in the T e g — P{ plane, for m = 1 . . .4 mixed modes of the stellar models introduced in Section |1] For each 
value of the azimuthal order, the extent of the re-mechanism instability at the rotation rates Q/Q c = 0.25 (dotted) and 0.5 (dashed) is 
indicated using thick lines. The small open triangles, filled squares, and open diamond show the observed periods of selected SPB stars, 
pulsating Be stars, and HD 208727, respectively (see text). Note that the ordinate scale is logarithmic. 



usually supposed (e.g., iRivinius et al-l EoOS. and references 
therein). This possibility is returned to in Section 

Having discussed the extent of the K-mechanism insta- 
bility, I now shift the focus to its strength. Fig. |5] plots the 
maximal growth rates of the instability across the HR di- 
agram, for the m = 1 ... 4 mixed modes at the rotation 
rate Q/Q c = 0.5. These values are intended to furnish some 
measure of the robustnes s of the instability. To this end, 
iDziembowski et al.l lll993fi employ the normalized growth 



rate introduced by IStellingwer i lll979l) . but in the absence 
of any sound theory of non-linear amplitude saturation, I 
prefer to adhere to the simple definition Q. 

The general trend seen in the figure is that the growth 
rates are at their greatest for the m = 1 modes, reaching 
up to a value rj ~ 1CF 4 ' 5 that is commensurate with the 
strongest instabilities associated with (3 Cep and SPB stars. 
Toward larger values of m, r\ declines rapidly, reflecting the 
fact that mixed modes having larger azimuthal orders ex- 
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hibit a more Rossby-like character (see, e.g. lLoul200Ct Fig. 1; 
note that the horizontal wavenumber k x plotted as abscissa 
is equivalent to m), and are therefore (Sec.|2l more difficult 
to excite via the k mechanism. 

A secondary trend in Fig. |K| is that the growth rates 
appear systematically larger for the stellar models situated 
near the TAMS. As shown in Fig. 4 of T05, these models 
are characterized by driving-zone thermal timescales that 
are rather longer than the corresponding ZAMS models. 
Presumably, mixed modes can satisfy the period matching 
requirement more easily at these long timescales, and there- 
fore are less susceptible to inhibition arising from radiative 
damping in the inner envelope. 

To bring the present section to a close, Fig. |S| plots the 
mixed-mode instability regions in the T e « — P plane. Here, 

Pi = -2tt/(ct - mfi) (8) 

is the pulsation period in an inertial frame of reference - a 
quantity of relevance to the observational community, since 
it can be ascertained directly from time-series monitoring of 
variable stars. The minus sign appears in this expression be- 
cause, for the mixed modes considered here, the co-rotating 
frame frequency a is always smaller than mfi. Also shown in 
the figure are the observed periods of the same stars plotted 
in Fig. [I] in cases where a star is multiperiodic, only the 
principal period is displayed. 

The figure reveals that the inertial frame periods be- 
come shorter with rotation rate, behaviour arising due to the 
appearance of in the denominator of the expression for 
Pi. More significantly, it is evident that the unstable m — 1 
mixed modes exhibit periods that are one or two orders of 
magnitude larger than those with azimuthal orders m 2. 
This is a consequence of the fact that the m = 1 modes have 
co-rotating frame frequencies that, to first order in Q, follow 
a ~ Q, (see eqn.0. Hence, the denominator in equation JHJ 
involves cancellation between two nearly equal terms, lead- 
ing to very large inertial- frame periods. 

As Fig. |S] demonstrates, these long periods are incom- 
patible with the much-shorter periods (~ 1 d) observed in 
SPB stars. While the mismatch between theory and observa- 
tions might be corrected by considering rotation rates above 
the present upper limit fl/Q c = 0.5, such rapid rotation 
would be difficult to reconcile with the angular frequencies 
£l/Sl c 5, 0-25 more typical to the majo rity of SPB stars (e.g., 
lAerts et al1ll999t lDeCat et alJl2005l) . The situation is not 
much improved for azimuthal orders m ^ 2, due to the fact 
that the instability regions at Q/S1 C = 0.25 - a more rep- 
resentative rotation rate - are situated at temperatures too 
cool to encompass many SPB stars. These findings amount 
to further evidence against a mixed-mode interpretation for 
the SPB phenomenon. 

The periods of the unstable m = 1 mixed modes in 
Fig.||j]are also much longer than those observed in pulsating 
Be stars. At the more-rapid rotation rates representative of 
these stars, it seems likely that this period discrepancy will 
be resolved. However, the m ^ 2 mixed modes, which match 
the typical observed periods quite well when Q/Q c = 0.5 
(even though the T c ff range of the instability is too cool), 
may then become inconsistent with the observations. This 
is one issue that will have to be addressed if mixed modes 
are considered a viable explanation for Be-star pulsation. 



6 DISCUSSION AND SUMMARY 

In the preceding sections, I use the non-adiabatic BOOJUM 
code to explore the stability of m = 1 ... 4 mixed modes 
in rotating B-type stars. The principal finding is that these 
modes are excited by the same iron-bump k mechanism re- 
sponsible for the pulsation of SPB and 8 Cep stars. This 
finding is significant, representing the discovery of a wholly- 
new type of instability in the early-type domain of the HR 
diagram. 

However, I should emphasize certain caveats relating 
to this discovery. As discussed in detail in T05, many ap- 
proximations are employed to derive the system of pulsation 
equations and boundary conditions solved by BOOJUM. Al- 
though these approximations are required for the problem 
to be tractable using present-day computational resources, 
they do mean that the results presented here should be re- 
garded as more qualitative than quantitative. 

In particular, neglect of the centrifugal deformation of 
the equilibrium stellar structure will certainly introduce a 
degree of error in the calculated eigenfrequencies. This is 
because r modes - and by extension mixed modes - tend 
to be rather sensitive to any departures from a sp herically- 
symmetric equilibrium state jProvost et al.lll98jJ) . Further- 
more, the results may be influenced by the incidence of dif- 
ferential rotation (re cently discovered in seismic analyses of 
two 8 Cep stars; see iDupret et aljl2004l iPamvatnvkh et aD 
2004), although it is as yet unclear what specific effect any 
departure from the assumed-uniform rotation might have. 

Nevertheless, while there certainly remains a great deal 
of scope for improving upon the theoretical foundation es- 
tablished in the present work, the case for the physicality 
of the mixed-mode instability in B-type stars is, I believe, 
convincing. Therefore, it is appropriate to discuss what im- 
plications it might have for studies of dynamical phenomena 
in early- type stars. 

As already made clear in the preceding section, I am 
skeptical that mixed modes can furnish a better explana- 
tion for the SPB phenomenon than the currently-accepted 
model of k mechanism excitation of (predominantly) I = 1 
g modes. The latter model has already proven very suc- 
cessful both in matching the o bserved distribution of SPB 
stars in the HR diagram (see iPamvatnvkbl ll999Tl . and in 
reproducing the variability exhibit ed by these stars (see 
lTownsendl|2002l : iDe Cat et alJl2005l) . Of course, this is not 
to say that there might not be some SPB stars that turn 
out to be mixed mode rather than g mode pulsators. In 
part icular, the so-call ed rapidly rotating Hipparcos SPBs 
(see lAerts et alJll999h may, upon closer scrutiny, be found 
to fall into this category. A thorough evaluation of such 
a possibility must await new tools for analyzing the spec- 
troscopic and photometric variations generated by mixed 
modes. These tools should prove relatively straightforward 
to derive from already-existi ng models based around the tra- 
ditional approximation (see iTowns end 1997. l2003bl) . 

The prospects for a mixed-mode role in pulsating Be 
stars seem rather more favourable. Near-critical rotation 
should shift the instability strips toward hotter temperatures 
and earlier spectral types, thereby allowing them to encom- 
pass a greater proportion of these stars. There may be diffi- 
culty in matching periods, and it has yet to be ascertained 
whether mixed modes can reproduce observed Be-star line 
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profile variations to the same degree that P o incare modes 
alrea dy have (see, e.g., iRivinius et alJ l200ll : iMaintz et alJ 
2003). However, in anticipation of the eventual resolution of 
these uncertainties, it is worthwhile reflecting on the part 
played by mixed modes in the episodic, disk-for ming mass 
loss from Be stars (see iPorter fc Riviniusll2003l and refer- 
ences t herein) . Following on f rom suggestions made bv lOsakil 
jl986h . lOwockil (12004 120051) has recently advanced a hy- 
drodynamical model for 'pulsationally driven orbital mass 
ejection' (PDOME), in which the addition of a pulsation- 
originated perturbation to the equatorial regions of a near- 
critical rotating star is sufficient to lift surface material into 
orbit. This material grad ually diffuses outward under the 
influence of viscosity (see iLee et al]ll99li : IPorte 3 ll99SL and 
references therein), to form a disk with the same Keplerian 
veloc ity distribution in f erred from obser vations of Be stars 
fe.g.. lDachs et alJll986tlHanuschii3ll99d) . 

As demonstrated bv lOwockil <l2005f) . a key ingredient in 
the success of the PDOME mechanism is that the azimuthal 
velocity perturbation is in the prograde direction when the 
density perturbation is maximal, ft would therefore seem 
that prograde pulsation modes are required for mass ejec- 
tion to occ ur, which is problem atical to reconcile with ob- 
servations dRivinius et al.ll2003t and references therein) that 
indicate retrograde-propagating modes in Be stars. However, 
an important distinction should be made here between group 
and phase velocities. While these velocities share the same 
sign in the case of Poincare, Kelvin and prograde mixed 
modes, the retrograde mixed modes behave differently: their 
phase velocity is (by definition) in the retrograde direction, 
yet their group velo city is oriented in the opposite direction 
(see, e.g., Gill 1982), with maximal density coinciding with 
prograde azimuthal velocity. Hence, they are able simulta- 
neously to fulfill both the observational constraint of retro- 
grade phase velocity, and the theoretical PDOME constraint 
of prograde group velocity. In tandem with a re-evaluation 
of t he degree to which Be stars approach critical rotation 
(see iTownsend et alJl2004) . retrograde mixed modes may 
therefore be the missing link required for Owocki's PDOME 
mechanism to operate. 

Looking now at the opposite, low-temperature end of 
the B spectral type, mixed modes may al so have a role 
to play in the phenomenon of Maia stars. IStruvd <ll955h 
first put forward the suggestion that there exists an in- 
stability strip in the HR diagram, characterized by pul- 
sation periods on the order of hours, that extends from 
late-B types to the early- A types 5 . Although the strip was 
named for its archetype Maia (20 Tau), subsequent ob- 
servations fe.g.. I&truve et all Il957t) of this particular star 
have, rather ironically, revealed it most likely to be non- 
variable. Other searches for stellar variabil i ty in the hypo- 
thetical Maia strip (e.g.. | McNamaralll985l : iLehmann et alJ 
ll995l:lPercv fc Wilsorj200cF ' have been largely unsuccessful. 
However, there are a handful of variable stars (e.g., 7 CrB - 
[Scholz et al1ll998t HD 208727 - iKallinger. Reeeen fc Weissl 
2002) that remain situated in the strip, and it would there- 



5 Man y authors quote the spectral range B7-A3, but IStruvd 
il955Tl was never this explicit, resorting instead to a sketched 
HR diagram. 



fore be premature to dismiss the Maia-star concept as mere 
will-o'-the-wisp. 

The discovery of the retrograde mixed-mode instability 
may provide a theoretical framework for understanding the 
origins of the elusive Maia stars. On the bas is of the stellar 
parameters quoted bv lKallinger et all i2002lh the Maia can- 
didate HD 208727 falls blueward of the ZAMS boundary in 
the HR diagram (see Fig.0J. However, if this boundary were 
shifted to hotter temperatures, for instance by considering a 
helium-enriched main sequence, the star would fall well in- 
side the m — 1 ... 4 mixed-mode instability regions. Further- 
more, the star's 7.6 h photometric period would be consis- 
tent with the inertial frame periods (Fig.|SJ typical to m = 3 
and m = 4 mixed modes for Q/Q c = 0.5 (in fact, the star 
probably rotates rat her faster than t his, fl /fl c > 0.67). Ac- 
cordingly, although Kall inger et alJ (120021) were dismissive 
of the ability of ordinary g modes to explain the variability 
of HD 208727, it appears that a mixed-mode interpretation 
may furnish a good match. 

The interpretation of Maia stars as mixed-mode pul- 
sators is lent ind i rect support in a recent paper by 
lAerts fc Kolenberd (l2005lh who point out that significant 
rotation - a prerequisite for the mixed-mode instability - 
appears to be a common factor in these stars. These authors 
then proceed to suggest that the Maia stars may simply be 
an extension of the SPB phenomenon to lower T e s arising 
due to the Coriolis force. The results presented in T05 would 
appear to rule out this latter possibility, by demonstrating 
(Fig. 3, ibid) that any Coriolis-originated extension to the 
SPB instability strip is insufficient to encompass the Maia 
stars. However, as Figs. [I] and [U of the present paper illus- 
trate, the mixed-mode instability at larger azimuthal orders 
(m > 4) can encompass cooler temperatures corresponding 
to early A spectral types, while at the same time exhibiting 
pulsation periods of a few hours. Although this establishes a 
case for a l ink between mixed modes and Maia stars, there 
should - as lAerts fc Kolenberd j2005l) themselves point out 
- be a significant improvement in observational datasets be- 
fore firm conclusions can be reached concerning this issue. 

Of course, irrespective of whether retrograde mixed 
modes have any connection with the Maia phenomenon, 
there are many pertinent questions that are raised by the 
discovery of their K-mechanism instability. What kinds of 
variability - spectroscopic or photometric - do these modes 
generate? Although the majority of SPB stars appear to be 
well understood as g mode pulsators, are there some that 
can be better better explained in terms of mixed mode exci- 
tation? Might mixed modes explain the variability and mass 
loss of pulsating Be stars? Can we expect mixed modes in 
other types of object, such as 7 Dor stars? I touch on some 
of these questions in the foregoing analysis, and hope further 
to address them, and others, in forthcoming papers. 



NOTE ADDED DURING REVISION 

A recent preprint bv lSavoniiei 1 2005i also investigates mixed- 
mode instability in early-type stars. Although there are some 
differences in the specific methodology and results, the gen- 
eral findings are in agreement with those reported in the 
present work. 
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